Nematics with quenched disorder : violation of self-averaging 



O 

(N 

(N 



43 

O 

B 

i 

-i— > 

» 



c 

o 
o 



> 

(N 

in 
q 

On 
O 

o 



X 



J. M. Fish and R. L. C. Vink 

Institute of Theoretical Physics, Georg-August-Universitat Gdttmgen, 
Friedrich-Hund-Platz 1, 37077 Gottingen, Germany 
(Dated: September 3, 2010) 

We consider the isotropic-to-nematic transition in liquid crystals confined to aerogel hosts, and 
assume that the aerogel acts as a random field. We generally find that self-averaging is violated. 
For a bulk transition that is weakly first-order, the violation of self-averaging is so severe, even the 
correlation length becomes non-self-averaging: no phase transition remains in this case. For a bulk 
transition that is more strongly first-order, the violation of self- averaging is milder, and a phase 
transition is observed. 
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Liquid crystals confined to quenched disordered me- 
dia are frequently encountered and of practical impor- 
tance In certain cases - the prototype example 
being silica aerogel - the disordered medium induces a 
quenched random field 2-8]. The random field couples 
to the liquid crystal at (essentially) arbitrary locations, 
and imposes a preferred orientation of the nematic direc- 
tor at these locations. One consequence of random field 
disorder in liquid crystals is the loss of long-range nematic 
order in all experimentally relevant dimensions d < 3 0- 
Il3j . This, however, does not rule out the existence of 
phase transitions. In contrast, the latter are routinely 
observed @, 0, H, EH and understanding the influence of 
random field disorder on liquid crystal phase transitions 
is an important topic. One known effect is that random 
fields can change the order of a transition 0, [l5[ . The 
bulk isotropic-to-nematic (IN) transition in three dimen- 
sions (3D) is usually first-order, but random fields can 
render this transition continuous 0, 0, [ljl [I?} or wipe 
it out completel y 3 - Other known effects include slow 
dynamics 19, Il8l420l| , lowering of phase transition tem- 
peratures [3J, |aL and formation of multidomain nematic 
structures [4|, [21] . 

It is also known that systems exposed to random fields 
generally do not self-average: results obtained for one 
sample of disorder, even if the sample is large, are not 
necessarily representative for all disorder samples [22i — 
|2~H . To what extent lack of self-averaging plays a role at 
the IN transition is the topic of the present Letter. Our 
main result is that, for a bulk IN transition in 3D that 
is weakly first-order, i.e. the experimentally most rele- 
vant case, the violation of self-averaging in the presence 
of random fields is so severe, even the correlation length 
becomes a non-self-averaging quantity [25|, [2|| . The IN 
transition temperature, as characterized by the temper- 
ature of the specific heat maximum, does not become 
sharp in the thermodynamic limit, but is given by a dis- 
tribution of finite width. Hence, no sharp phase transi- 
tion remains. 

To illustrate this point, we have simulated the sprin- 
kled silica spin (SSS) model [27|; models such as this are 
routinely used to describe nematics in disordered media 
OS [U M^k- The SSS model is defined on a 3D peri- 



odic V = L x L x L lattice. A 3D unit vector di (spin) 
is attached to each lattice site i. The energy density is 
given by 



-J/Vj2\di-dj\ p , J>0, 



(1) 



with the sum over nearest neighbors (in what follows, 
the temperature T is expressed in units of J/ks, with 
ks the Boltzmann constant). We set p = 2 for now; 
Eq.([T]) then resembles the Lebwohl-Lasher model [30| . 
which undergoes a weak first-order IN transition from 
a high-T isotropic phase (exponential decay of the ne- 
matic correlation function to zero), to a low-T nematic 
phase with long-range order (exponential decay of the ne- 
matic correlation function to a finite positive value). In 
the SSS model, quenched disorder is introduced by mark- 
ing a fraction q of randomly selected spins as quenched 
(we use q = 0.1 always). These spins are oriented ran- 
domly at the start of the simulation and remain static 
thereafter, which can be conceived as a random field of 
infinite strength acting on a fraction of the spins. Even 
though the random field strength is infinite, q = 0.1 re- 
mains in the weak field limit, in the sense that the non- 
quenched spins still form a percolating network. If q is 
set above the percolation threshold, any phase transition 
gets trivially blocked, since then the correlations cannot 
propagate through the lattice anymore. The SSS model 
is different from the random-field Ising model because the 
spins are 3D continuous vectors, as opposed to discrete 
integers. The SSS model does not support long-range 
nematic order at any finite temperature [9(. 

Most of our analysis is based on the distribution 
Pj^fL{e, s), defined as the probability to observe energy 
density e and nematic order parameter s, at tempera- 
ture T, system size L, and for some sample of random 
fields k. We measure the distributions for L = 7 — 15. 
The nematic order parameter s is defined as the maxi- 
mum eigenvalue of the 3D orientational tensor. In a per- 
fectly aligned nematic sample s = 1, while an isotropic 
sample yields s — > in the thermodynamic limit. We use 
broad histogram methods, namely Wang-Landau sam- 
pling [3l| and successive umbrella sampling (32j, to ob- 
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FIG. 1: Rs versus T using p = 2 (a) and p — 10 (b) for several 
L. The temperature where R s is maximal defines Tr. Note 
that Tr decreases with L. For p — 2, there is no self-averaging 
at low T. In contrast, for p = 10, self-averaging is restored at 
low T, and a sharp phase transition occurs (marked with the 
dot). 




0.10 1 1 1 1 

5 7 10 15 

r 



(b) 1.10 




5 7 10 15 



FIG. 2: Correlation functions for p = 2, L = 30, and T = 0.5 
(which is well below T/j) on double logarithmic scales; due 
to periodic boundaries up to r max = 15 can be sampled, (a) 
G(r) obtained for several samples (dashed curves) together 
with the disorder-averaged result [G(r)\ (solid curve), (b) k 
versus r; the dashed line is a power law fit to the large r 
regime. 



tain P£ j,(e, s). These methods ensure that the simula- 
tion performs a random walk in phase space. This is 
crucial because the SSS model is known to exhibit meta- 
stable states [2^, [3^] , in which standard Monte Carlo sim- 
ulations (sampling directly from the Boltzmann distribu- 
tion) may "get stuck" . Since we expect self-averaging to 
be violated, it is crucial that the distributions be mea- 
sured for k = 1, . . . , M random field samples, where M 
must be large. We use M ~ 1000 — 2500, based on the 
convergence of "running averages" of quantities of in- 
terest onto plateau values. We also measure correlation 
functions; the latter are obtained for L = 30 using stan- 
dard Boltzmann sampling. 

For each sample k, we compute the thermally averaged 
nematic order parameter (s}k and measure the fluctua- 
tion between samples R 2 = [(s) 2 ] — [(s)] 2 , with [•] the 
disorder average [X n ] = (1/M) YlkLi If the system 
self-averages, R s — > in the thermodynamic limit, in 
which case a single experiment on a large system will be 
representative for all samples. In Fig. [Ha), we plot R s 
versus T for three system sizes. The striking result is 
that, at low temperature, R s does not decay to zero with 
increasing L but remains finite. The onset to the non-self- 
averaging regime is marked by a maximum in i? s , at tem- 
perature T = Tr. We thus identify two regimes: a high-T 
regime (T > Tr) where the SSS model self-averages (i? s 
decreases with L), and a low-T regime (T < Tr) where 
self- averaging is violated (R s remains finite). 

The violation of self-averaging at low T profoundly af- 
fects the nematic correlation function G(r) — (|(d(0) ■ 

d(r)) 2 - i) E3- (In this work, G(r) is calculated using 
all spins, i.e. free and static ones.) Since it holds that 
G(L/2) — (s) 2 , with L the edge of the simulation box, 



and since R s > 0, fluctuations in G(r) between disor- 
der samples are automatically implied. We must there- 
fore consider Gk(r), i.e. the nematic correlation function 
obtained in the fc-th random field sample. In the high- 
T regime, we find that Gfc(r) decays exponentially to 
zero, with negligible fluctuations between samples: the 
SSS model is isotropic and self-averaging when T > Tr. 
In contrast, in the low-T regime, Gfe(r) fluctuates pro- 
foundly between disorder samples (Fig. [Ha)). Note that 
we concentrate on the tail of G{r) and so the range r < 5 
is discarded. In some samples, Gk(r) decays very rapidly, 
while in others the decay is much slower. Clearly, when 
T < Tr a single measurement of Gk(r) is not represen- 
tative. 

The key point is that, in random field systems, there 
exist two correlation functions: the connected correlation 
function [G(r)] (i.e. the nematic correlation function av- 
eraged over all samples), and the disconnected correlation 
function [G(r) 2 } @ [35|. The solid curve in Fig. [2fa) 
shows [G(r)}: its decay to zero is most consistent with 
a power law, suggesting quasi-long-range order. This 
agrees with Ref. [l3|, but it disagrees with Ref. [13 (where 
short-ranged exponential decay is observed). Regardless 
of the precise form of the decay, we confirm that Gk(r) 
does not self-average. This is shown in Fig. [2jb), where 
k = [G(r) 2 ]/[G(r)] 2 is plotted. At large r, power law 
growth k oc t® , with 9 ~ 0.1, is observed. The discon- 
nected correlations thus decay independently from the 
connected ones. In contrast, if Gfc(r) were self-averaging, 
the fluctuation [G(r) 2 ] — [G(r)\ 2 would be zero at large 
r: [G(r) 2 } and [G(r)] 2 then decay with the same expo- 
nent. Since the correlation functions do not self-average, 
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FIG. 3: Histograms of T c ^ (a) and c maXi fc (b), shifted by their 
respective averages, and for several L. The histograms do not 
become sharp as L increases. 



it follows that properties extracted from these functions 
do not self-average either, which includes the correlation 
length £ [IE, EH . The Brout argument (3(| , which con- 
ceives the thermodynamic limit as a large number of in- 
dependent sub-samples of size £, thus breaks down. In- 
stead, £ must be regarded as a random variable. The 
power law decay of [G(r)] observed by us indicates that 
£ itself is very large, if not infinite. For £ — > oo, the Brout 
argument breaks down in any case [24j . 

How does this affect the IN transition in the SSS 
model? The usual approach to detect the IN transi- 
tion is to measure the specific heat c = V((e 2 ) — (e) 2 ) 
versus T; at the transition, c reaches a maximum. For 
each random field sample k, we measured the temper- 
ature T C: fc where c was maximal, and the correspond- 
ing value c maXj fe. Since £ does not self-average, an un- 
usually large fluctuation [T 2 ] — [T c ] 2 is expected. This 
is confirmed in Fig. [SJa), where histograms of T Ct u are 
shown, shifted such that [T c ] is at zero, and for several 
L. The striking result is that the distributions do not 
become sharp as L increases. The specific heat itself is 
also non-self-averaging. This is illustrated in Fig. Elb), 
where histograms of c maXj fc are shown, shifted by [c max ], 
and again for several L. We also observed that [T c ] is 
very close to the temperature Tr where R s is maximal. 
A signature of the onset to the low-T regime (where self- 
averaging is violated) is thus also provided by the specific 
heat maximum. Both [T c ] and Tr decrease with increas- 
ing L: the non-self-averaging regime T < Tr thus gets 
smaller in larger systems. Unfortunately, finite size scal- 
ing with a non-self-averaging correlation length is compli- 
cated - a rigorous scaling theory remains elusive - and so 
it is difficult to estimate Tr in the thermodynamic limit. 
The decrease of Tr with L, and hence of [T c ], is in any 
case slow. For instance, if we assume a power law shift 



Tr — Too oc 1/L y , Too = liniL^oo Tr, a fit to our data 
yields a maximum value for the exponent !/ max ~ 0.16; 
this upper bound is obtained by assuming = 0. 

To conclude: the SSS model (with p = 2 in Eq.Q 
and quenched spin fraction q = 0.1) does not feature a 
sharp phase transition. For a given sample k of random 
fields, a well-defined temperature T c ^ where the specific 
heat attains its maximum can be measured, but the fluc- 
tuation in T c ^ between samples remains finite, even as 
L — > oo. We attribute this behavior to the existence of a 
non-trivial disconnected correlation function, which im- 
plies a non-self-averaging correlation length when T < 
Tr ~ [T c ]. In this regime, the SSS model does not self- 
average. The temperature Tr decays extremely slowly 
with system size; whether Tr remains finite in the ther- 
modynamic limit, or whether it decays to zero, cannot 
be discerned from our data. Since the decay of Tr and 
[T c ] with L is slow, it is likely that the non-self-averaging 
regime survives in macroscopic samples (even if Too = 0) . 
We expect that by varying T a maximum in the specific 
heat will be found, but the value of the specific heat at 
the maximum will vary between samples. There is some 
experimental evidence for this behavior. The liquid crys- 
tal 8CB in bulk undergoes a weak first-order IN transi- 
tion 0, Q, as does Eq.(p} with p — 2. Upon insertion in 
aerogel, the enthalpy obtained in different samples ranges 
from 3.6 — 5.23 J/g, which is unusually large 0j- How- 
ever, since the enthalpy is related to the specific heat, and 
since the specific heat does not self-average (Fig. [3^b)), 
a large enthalpy fluctuation between samples would, in 
fact, not be unexpected. 

Do our results imply the absence of IN transitions, in 
general, in the presence of random field disorder? The 
answer to this question is an unequivocal "No" ! The 
phase behavior of liquid crystals is not dictated by any 
universality class, and by changing details in the parti- 
cle interaction qualitatively different scenarios may de- 
velop To illustrate this, we reconsider Eq.fll]) using 
p = 10; this makes the bulk IN transition more strongly 
first-order (37j . Again using a fraction of quenched spins 
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q = 0.1, we show in Fig. QJb) the variation of R s with T. 
The striking difference with p — 2 is that self-averaging is 
restored at low temperatures. There now appears an in- 
termediate regime of temperatures where self averaging is 
violated, but this regime becomes smaller as L increases. 
Hence, in the thermodynamic limit, self-averaging is vio- 
lated at only one temperature, which then reflects a sharp 
phase transition, with finite size effects given by 22] 

VP?} - [T c ] 2 oc [T e ] - ^ oc 1/L» (2) 

where Too is the transition temperature in the thermo- 
dynamic limit. This implies that histograms of (Tx> — 
T c ,k)L v become T-independent, provided correct values 
of Too and y are used. The scaling is confirmed in Fig. 21 
using Too ~ 0.558 and y m 0.88, and the collapse is clearly 
excellent. Incidentally, Tx) corresponds to an approxi- 
mate intersection in curves of R s versus T for different 
L (Fig. Hfb)), which offers an alternative route to locate 
the transition. 

In summary: we have shown that the IN transition 



in the presence of random fields is strongly affected by 
a lack of self-averaging. Certainly for computer simula- 
tions, taking a disorder average [•] involving many sam- 
ples is crucial. For a bulk IN transition that is weakly 
first-order, the violation of self-averaging is so severe, 
even the correlation length £ becomes non-self-averaging 
[25L I2H ] . This manifests itself from the nematic correla- 
tion function, which becomes strongly sample dependent. 
A consequence is that no sharp IN transition remains in 
this case. For a bulk IN transition that is more strongly 
first-order, the violation of self- averaging is restricted to 
a single temperature in the thermodynamic limit. In this 
case, a phase transition does occur, and finite size effects 
near the transition are well understood 22]. As far as 
we know, a scaling theory for the case where £ does not 
self-average remains elusive. In some sense, the finite size 
effects observed by us for p = 2 resemble those of Eq. ([2]) , 
but in the limit where y — > 0. Perhaps a new scaling 
theory should be developed keeping this in mind. 

This work was supported by the Deutsche Forschungsge- 
meinschaft, Emmy Noether program: VI 483/1-1. 



[1] Crawford GP and Zumer S (Eds.), Liquid Crystals in 
Complex Geometries (Taylor & Francis, 1996) 

[2] Wu L, Zhou B, Garland CW, Bellini T, and Schae- 
fer DW, Heat-capacity study of nematic-isotropic and 
nematic- smectic- A transitions for octylcyanobiphenyl in 
silica aerogels, Phys. Rev. E 51, 2157 (1995) 

[3] Maritan A, Cieplak M, and Banavar JR, Nematic- 
Isotropic Transition in Porous Media, in: GP Crawford 
and S Zumer (Eds.), Liquid Crystals in Complex Geome- 
tries, 483 (Taylor & Francis, 1996) 

[4] Bellini T, Clark NA, Muzny CD, Wu L, Garland CW, 
Schaefer DW, and Oliver BJ, Phase behavior of the liquid 
crystal 8CB in a silica aerogel, Phys. Rev. Lett. 69, 788 
(1992) 

[5] Wu XI, Goldburg WI, Liu MX, and Xue JZ, Slow Dy- 
namics of Isotropic-Nematic Phase Transition in Silica 
Gels, Phys. Rev. Lett. 69, 470 (1992) 

[6] Iannacchione G and Finotello D, Confinement and 
orientational study at liquid crystal phase transitions, 
Liq. Cryst. 14, 1135 (1993) 

[7] Tripathi S, Rosenblatt C, and Aliev FM, Orienta- 
tional susceptibility in porous glass near a bulk nematic- 
isotropic phase transition, Phys. Rev. Lett. 72, 2725 
(1994) 

[8] Zhou B, Iannacchione GS, Garland CW, and Bellini 
T, Random-field effects on the nematic-smectic-A phase 
transition due to silica aerosil particles, Phys. Rev. E 55, 
2962 (1997) 

[9] Imry Y and Ma SK, Random-Field Instability of the Or- 
dered State of Continuous Symmetry, Phys. Rev. Lett. 
35, 1399 (1975) 

[10] Bellini T, Buscaglia M, Chiccoli C, Mantegazza F, Pasini 
P, and Zannoni C, Nematics with Quenched Disorder: 
What Is Left when Long Range Order Is Disrupted?, 
Phys. Rev. Lett. 85, 1008 (2000) 

[11] Cleaver DJ, Kralj S, Sluckin TJ, and Allen MP, The Ran- 
dom Anisotropy Nematic Spin Model, in: GP Crawford 



and S Zumer (Eds.), Liquid Crystals in Complex Geome- 
tries, 467 (Taylor & Francis, 1996) 

[12] Feldman DE and Pelcovits RA, Liquid crystals in ran- 
dom porous media: Disorder is stronger in low-density 
aerosils, Phys. Rev. E 70, 040702 (2004) 

[13] Feldman DE, Quasi-Long-Range Order in Nematics Con- 
fined in Random Porous Media, Phys. Rev. Lett. 84, 4886 
(2000) 

[14] Clark NA, Bellini T, Malzbender RM, Thomas BN, Rap- 
paport AG, Muzny CD, Schaefer DW, and Hrubesh L, X- 
ray scattering study of smectic ordering in a silica aerogel, 
Phys. Rev. Lett. 71, 3505 (1993) 

[15] Imry Y and Wortis M, Influence of quenched impurities 
on first-order phase transitions, Phys. Rev. B 19, 3580 
(1979) 

[16] Iannacchione GS, Crawford GP, Zumer S, Doane JW, 
and Finotello D, Randomly constrained orientational or- 
der in porous glass, Phys. Rev. Lett. 71, 2595 (1993) 

[17] Kralj S, Lahajnar G, Zidansek A, Kopac NV, Vilfan M, 
Blinc R, and Kosec M, Deuterium NMR of a pentyl- 
cyanobiphenyl liquid crystal confined in a silica aerogel 
matrix, Phys. Rev. E 48, 340 (1993) 

[18] Mercuri F, Paoloni S, Zammit U, and Marinelli M, Dy- 
namics at the Nematic-isotropic Phase Transition in 
Aerosil Dispersed Liquid Crystal, Phys. Rev. Lett. 94, 
247801 (2005) 

[19] Bellini T, Clark NA, and Schaefer DW, Dynamic Light 
Scattering Study of Nematic and Smectic- A Liquid Crys- 
tal Ordering in Silica Aerogel, Phys. Rev. Lett. 74, 2740 
(1995) 

[20] Rotunno M, Buscaglia M, Chiccoli C, Mantegazza F, 
Pasini P, Bellini T, and Zannoni C, Nematics with 
Quenched Disorder: Pinning out the Origin of Memory, 
Phys. Rev. Lett. 94, 097802 (2005) 

[21] Petridis L and Terentjev EM, Nematic-isotropic transi- 
tion with quenched disorder, Phys. Rev. E 74, 051707 
(2006) 



■5 



[22] Aharony A and Harris AB, Absence of Self-Averaging and 

Universal Fluctuations in Random Systems near Critical 

Points, Phys. Rev. Lett. 77, 3700 (1996) 
[23] Wiseman S and Domany E, Finite-Size Scaling and 

Lack of Self-Averaging in Critical Disordered Systems, 

Phys. Rev. Lett. 81, 22 (1998) 
[24] Wiseman S and Domany E, Self-averaging, distribution 

of pseudocritical temperatures, and finite size scaling in 

critical disordered systems, Phys. Rev. E 58, 2938 (1998) 
[25] Parisi G, Picco M, and Sourlas N, Scale invariance and 

self- averaging in disordered systems, Europhys. Lett. 66, 

465 (2004) 

[26] Parisi G and Sourlas N, Scale Invariance in Disordered 
Systems: The Example of the Random-Field Ising Model, 
Phys. Rev. Lett. 89, 257204 (2002) 

[27] Bellini T, Chiccoli C, Pasini P, and Zannoni C, 
Lattice Spin Models of Liquid Crystals in Aerogels, 
Mol. Cryst. Liq. Cryst. 290, 227 (1996) 

[28] Chakrabarti J, Simulation Evidence of Critical Behav- 
ior of Isotropic-Nematic Phase Transition in a Porous 
Medium, Phys. Rev. Lett. 81, 385 (1998) 

[29] Maritan A, Cieplak M, Bellini T, and Banavar 
JR, Nematic-Isotropic Transition in Porous Media, 
Phys. Rev. Lett. 72, 4113 (1994) 



[30] Lebwohl PA and Lasher G, Nematic-Liquid-Crystal Or- 
der: A Monte Carlo Calculation, Phys. Rev. A 6, 426 
(1972) 

[31] Wang F and Landau DP, Efficient, Multiple-Range Ran- 
dom Walk Algorithm to Calculate the Density of States, 
Phys. Rev. Lett. 86, 2050 (2001) 

[32] Virnau P and Miiller M, Calculation of free energy 
through successive umbrella sampling, J. Chem. Phys. 
120, 10925 (2004) 

[33] Buscaglia M, Bellini T, Chiccoli C, Mantegazza F, Pasini 
P, Rotunno M, and Zannoni C, Memory effects in ne- 
matics with quenched disorder, Phys. Rev. E 74, 011706 
(2006) 

[34] Nattermann T, Theory of the Random Field Ising Model, 
in: AP Young (Ed.), Spin Glasses and Random Fields, 
277 (World Scientific, Singapore, 1998) 

[35] Khasanov B, Isotropic phase of nematics in porous media, 
JETP Lett. 81, 24 (2005) 

[36] Brout R, Statistical Mechanical Theory of a Random Fer- 
romagnetic System, Phys. Rev. 115, 824 (1959) 

[37] Fish JM and Vink RLC, Isotropic-to-nematic transition 
in confined liquid crystals: An essentially nonuniversal 
phenomenon, Phys. Rev. E 81, 021705 (2010) 



